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ABSTRACT 


A test was conducted to determine the longitudinal dynamic characteristics of the 
Stub D Atlas without propellants or payload. The dynamic characteristics obtained are 
given in terms of modal damping ratios, vehicle mode shapes, and acceleration responses 
to a force input over the frequency range of 1 to 120 Hz, Natural frequencies were deter- 
mined by comparing the acceleration response results with the mode shape and ortho- 
gonality test results. The results showed that a one-dimensional model may not be suf- 
ficient to predict the behavior of the vehicle. 
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SUMMARY 

A test was conducted to determine the longitudinal dynamic characteristics of the 
Stub D Atlas without propellants or payload. The dynamic characteristics obtained are 
given in terms of modal dampii^ ratios, vehicle mode shapes, and acceleration responses 
to a force input over the frequency range of 1 to 120 hertz. Natural frequencies were 
determined by comparii^ the acceleration response results with the mode shape and ortho- 
gonality test results. The results showed that a one-dimensional model may not be suf- 
ficient to predict the behavior of the vehicle. 


INTRODUCTION 

A definition of launch vehicle longitudinal dynamics is required in the study of vibra- 
tions caused by dynamic coupling between structure, ei^ines, and propellant lines (com- 
monly referred to as POGO) and to determine vehicle structural loads. The dynamic 
characteristics desired are in terms of modal damping ratios, vehicle mode shapes, and 
acceleration responses at various vehicle locations to a force input. 

An experimental study of the longitudinal dynamics of a Stub D Atlas with water as a 
simulated propellant, together with the Centaur vehicle and the Surveyor spacecraft as 
payload, has been reported by Gerus and Housley in reference 1. They also studied the 
liquid-oxygen -duct inlet pressure response. The experimental results obtained were 
compared with General Dynamics/Convair's theoretical values. These theoretical values 
were based on a one-dimensional lumped spring-mass model of the system which treats 
fluids in the manner described in reference 2. The analysis based on this model pre- 
dicted natural frequencies reasonably well but did not predict the amplitude of the pres- 
sure oscillations at the liquid-oxygen-duct inlet (tank base) with adequate accuracy. This 
pressure is one of the key factors in the analysis of the PCX50 oscillations. 



From the results of reference 1, it was evident that an improvement of the present 
one -dimensional model of the vehicle was in order. Before this problem can be attacked, 
a better description of the empty Atlas djTiamic response and its contribution to the over- 
all damping is needed (ref. 3). This investigation was conducted to determine the longi- 
tudinal dynamic characteristics of a Stub D Atlas vehicle without payload or propellants. 


TEST APPARATUS AND PROCEDURE 
Test Vehicle 

The test vehicle employed was a Stub D Atlas S/N 116 with a Centaur interstage 
adapter attached, A stiffening rii^ was added to the top of the interstage adapter to allow 
attachment of lateral stabilization springs. 

The booster and sustainer ermines erf the vehicle were removed and replaced with an 
equivalent-weight X-frame which served also as anchor points for the suspension system 
(fig. 1). The X-frame mass was 5400 pounds mass (2449 kg) and was considered to be 
rigid over the frequency range of interest. 

Durir^ all tests, the Atlas was maintained at standby pressures of 8 and 16 psig 
(5. 516x10^ and 11. 032x10^ N/m^ gage) for the liquid-oxygen and fuel tanks, respectively. 
The tanks contained no liquids. These pressures were chosen so that personnel could 
enter the testing area durii^ testii^ without risk. 


Suspension System 

The vehicle was supported during the dynamic tests by a suspension system (fig. 1) 
which consisted of a rigid frame at the base of the vehicle (X-frame) suspended by four 
steel cables, each of which was in him suspended by a sprir^ box containii^ 16 springs 
each (fig. 2). The sprii^ boxes were, in turn, connected to the tower strueture through 
a liftir^ and vehicle -weight-measuring system. 

Lateral stability was provided by a set of four radial springs attached to the stiffen- 
ing rirg at the top of the interstage adapter and a set of eight tangential springs on the 
X-frame, The spring rates and arrangements are shown in figures 3 and 4. 

Force input to the Atlas was effected at the booster engine attachment points on the 
vehicle y-axis. The vehicle also contacted the X-frame at the x-axis; however, these 
points served only to stabilize the vehicle and, hence, carried only a small part of the 
total weight. 

The suspension system spring constant was chosen so that the rigid-body mode 
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would be within the measurable frequency range. The total stiffness of the suspension 
system was set (by using 16 springs per box) at 25 866 pounds per inch (4. 530x10® N/m). 
The effective stiffness, however, was calculated from the observed natural frequency at 
20 500 pounds per inch (3. 590x10 N/m). This was so because a number of individual 
sprii^s were not always engaged due to variations in the unloaded spring lengths. The 
exciter flexures had a stiffness (in the axial direction) of 7000 pounds per inch (1.226x10® 
N/m), which becomes part of the suspension system stiffness during logarithmic decay. 

Instu mentation 

Three basic types of measurements were made: skin accelerations, force input from 
%he exciter, and displacements. 

Piezoelectric -type accelerometers were used for all but the low-frequency (f < 10 Hz) 
^testing. The accelerometers were mounted on small plastic cubes which were in turn 
glued to the vehicle skin (fig. 5). All accelerometers were mounted on the skin of the 
vehicle so as to detect longitudmal motion. Twenty-eight accelerometers were mounted 
on the +x-axis side and six on the +y-axis side of the vehicle, as shown in figure 6. 

Exciter force input to the X-frame was measured with a strain-type 15 000-pound 
(6. 672xlO^-N) load cell mounted between the exciter and the X-frame (see fig. 7). The 
piezoelectric-type accelerometer and load cell measurement accuracies were estimated 
to be better then 4 percent and 1 percent of full scale, respectively. 

For the low-frequency testing (f < 10 Hz) of the vehicle suspension system, linear 
potentiometers were used to indicate X-frame and interstage adapter position. At these 
frequencies, the linear potentiometers exhibited no dynamic degradation. Also, ten 
strain-gage-type accelerometers, with flat low-frequency response, were distributed 
over the vehicle, as shown in figure 8. 

The vehicle was driven aloi^ the longitudinal axis by an electrodynamic shaker 
through the load cell and the X-frame (fig. 9). The suspension system was used to 
covinteract the weight of the vehicle and, thus, aUow the shaker to supply only the dynamic 
part of the force . For high frequency ^ > 10 Hz), the shaker current was supplied by an 
electronic amplifier. The test procedure for high frequency was simply to drive the 
system at the frequency of interest until sinusoidal steady-state was achieved and record 
a suitable sample of the instrument response through a digitizer system. For all sinu- 
soidal tests, the driving force was adjusted so that the acceleration never exceeded 1 g 
(9. 807 m/sec^) on the vehicle or ±l/4-inch (0. 635xl0"^-m) displacement of the X-frame. 

o 

GeneraUy, acceleration amplitude was kept at the 3/4-g (7. 355-m/sec ) level. Under 
these conditions, the force levels ranged from 2000 to 13 000 pouiwis (8. 896x10® to 
57. 827x10® N). 

The test procedure to determine the vehicle suspension system characteristics at low 
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frequency varied from that at high frequency in that an amplidyne was used to supply 
exciter current. 

To find the dampii^ of a given mode, the system was driven at the modal frequency 
until sinusoidal steady-state was achieved. The shaker power was then cut off and the 
decay of a representative set of vehicle accelerations at 12 vehicle positions was record- 
ed on magnetic tape (analog signals). The dampii^ ratios were then determined by using 
the logarithmic decay technique. 


TEST RESULTS AND DISCUSSION 

The magnitude and phase angle of the acceleration relative to force of 28 positions on* 
the 4 -x-axis side and 6 positions on the +y-axis side of the vehicle are plotted against 
frequency and presented in figures 10 and 11, ^ 

These responses were obtained from the steady -state sinusoidal tests. The fre- 
quency range was from 1 to 120 hertz. The acceleration-relative -to-f or ce plots presented 
herein show the response from 10 to 120 hertz. These plots show a great number of 
resonances, but the most predominant are at 19, 40, 47, 72, 78, and 117 hertz. The 
response of the displacement relative to force of the X-frame is plotted against fre- 
quency in the 1- to 12. 5-hertz range (fig, 12). The suspension system vehicle mode 
(rigid-body mode) for the system occurred at 3. 5 hertz. 

Other than the rigid-body natural frequency, the determination of the natural fre- 
quencies from the plots of amplitude ratio against frequency proved to be a rather dif- 
ficult task because of the number of resonances that were observed. For this reason, 
vehicle djmamic shape plots of acceleration magnitude and phase angle relative to force, 
both against station, were obtained for a number of frequencies (fig. 13). These plots 
show that the loi^itudinal responses of all stations were nearly in phase at the lower fre- 
quencies (relative to 47 Hz), As the frequencies increased, the phase gradually shifted 
until, at 47 hertz, there was a 180° phase shift from the base of the vehicle to the top. 

This indicates that the resonances v4iich occurred below 47 hertz were oriented in other 
planes and may have been coupled, lateral, or shell resonances. From these consider- 
ations, it may be concluded that 47 hertz is the first vehicle longitudinal mode. 

Orthogonality tests were conducted for a number of resonance points and the results 
are presented in table I. The normalized orthogonality condition is given by 
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where l is the length of the vehicle, W(z) is the mass per xinit length of the vehicle, and 
Uj^(z)/F and ’u^(z)/F are the dynamic shapes of the vehicle at resonant test frequencies 
n and m. This tjq)e of orthogonality normalization is given in reference 4. The inte- 
grals were numerically evaluated using Gauss' quadrature method. 

In order to apply the orthcsonalilw condition, the following assumptions are made: 

(1) The effect of vehicle damping was negligible. 

(2) Since the first accelerometer was at station 548 and not at the top of the vehicle 

« (station 413), the weight of the vehicle from stations 413 to 548 was considered to be part 
of the weight of the first summation increment, 

(3) The X-frame was assumed to be rigid. Therefore, the bottom of the X-frame 
(station 1278) has the same dynamic response as the top (station 1254). 

The orthogonality expression was evaluated by using the digital computer. The 
least-squares method was used in curve fitting the dynamic shapes from known data 
points. A plot of vehicle mass distribution against distance is presented in figure 14. 

The near orthogonality between the 3, 5- and 47 -hertz modes indicates that 47 hertz 
is the first vehicle longitudinal mode. This agrees with the mode-shape and phase -angle- 
against- station results. The results also showed tiiat tiie 92. 8-hertz mode was ortho- 
gonal to the 47 -hertz mode but was not orthogonal to the 3, 5-hertz mode. 

The amplitude ratio plots showed rather large resonance at 72, 78, and 117 hertz. 
The vehicle dynamic plots show that the phase angle at 72 and 78 hertz shifts through 
285° and 305°, respectively, in comparison to 360° for a second flexible mode, while it 
is inconclusive at 117 hertz. The phase angle at 110 hertz, however, shifts through 
540°, which is characteristic of the third flexible mode, Indicating that the 117-hertz 
resonance is beyond the range of the second mode. The orthogonality test results did not 
show either the 72- or the 78-hertz resonance as orthogonal to the -3. 5- and 47-hertz 
modes . 

The y-axis side results for stations 572, 596, 764, and 932 (figs. 11(a) to (d)) match 
their respective x-axis results (figs. 10(b), (c), (e), and (h)) in both amplitude and phase 
angle. The y-axis side results for stations 1148 and 1176 (figs. 11(e) and (f)) do not 
match their respective x-axis side results (figs. 10(n) and (o)) at all frequencies. 

For stations 1148 and 1176, the amplitude ratio on the y-axis side at 47. 5 hertz was 
larger than the x-axis side amplitude ratio by factors of 2. 80 and 2. 20, respectively. At 
higher frequencies, the reverse was true; x-axis values exceeded y-axis values at 78, 83, 
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and 92. 5 hertz by an average factor of 1.95, 1. 74, and 1.21, respectively. The phase 
angle was different between the two sides from 36 to 70 hertz and from 100 to 120 hertz. 
This difference indicates that the dynamic characteristics are the same for the upper 
portion of the vehicle for both x- and y-axis sides but are not the same for the lower por- 
tion. This can be attributed to the fact that the upper portion of the vehicle is symmet- 
rical while the lower portion is not. The asymmetry is caused by the long and short pods 
attached to the skin of the vehicle on the -y and +y sides, respectively. 

Damping was calculated by using the logarithmic decay method, the dampii^ ratio 
^ given by 


^ = 


1 

JLln_£ 
2 itN a 


N 


where ? is the ratio of the actaal damping to the critical dampii^ and A- and A,, are 

til ^ ^ 

the amplitudes at the zero and N cycles of the decaying sinusoidal wave. The damping 

associated with the suspension system was 5 percent of critical based on Ic^arithmic 
decay and 4. 6 percent of critical based on the peak value of the amplirude. The damping 
values at 47 and 72 hertz were 1. 6 and 1, 45 percent (rf critical, respectively, based on 
Ic^arithmic decay. These modal damping constants are functions of the vehicle and sus- 
pension system local damping. The suspension system dampii^, in particidar, can 
greatly influence the first flexible mode damping even though the frequencies are well 
separated. The relative magnitude of this effect can be determined by the methods of 
reference 3. Some of the modal frequencies, percentages of critical damping, and their 
respective transient decay plots are shown in figure 15. 


CONCLUDING REMARKS 

From the study of the orthogonality test results, vehicle dynamic shapes, and fre- 
quency responses, we can reasonably conclude that the rigid-body mode of the vehicle and 
suspension system is at 3. 5 hertz and the first vehicle longitudinal flc'xible mode is at 
47 hertz. From the amplitude ratio results and the vehicle dynamic shapes the 78-hertz, 
or possibly the 72 -hertz, mode might be considered to be the second flexible mode. How- 
ever, this possibility could not be substantiated by orthogonality results. Damping at 
these modes, using the logarithmic decay method, was 5.0, 1.6, and 1.45 percent of 
critical, respectively. 

Since results for the x-axis and y-axis sides compare extremely well for stations 
572 to 932, it can be concluded that the d 3 mamic characteristics of the upper portion of 
the vehicle are the same for both sides over the test frequency range. 
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For the lower portion of the vehicle, the x- and y-axis results did not match, which 
means that the dynamic characteristics there are not the same for both sides. This is so 
because the vehicle is asymmetrical at the lower portion. In the upper part of the vehicle 
where the results matched, the vehicle is symmetrical. 

Because of the number of resonances shown in the acceleration response results, a 
one-dimensional spring-mass model may not be sufficient to predict the behavior of the 
vehicle. The apparent coupling of modes and the number of resonances indicates that 
more refined models incorporating bendii^ and shell effects might be necessary for empty 
and near- empty conditions , 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 9, 1968, 

180-31-01-01-22. 
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APPENDIX - SYMBOLS 


A amplitude of sinusoidal wave 

F applied force, lb (N) 

f frequency, Hz 

I length of vehicle, in. (m) 

u longitudinal displacement, in. (m) 

u longitudinal acceleration, g 

(m/sec^) 

W mass per unit length of the vehicle 
mass/in. (kg/m) 


X horizontal direction as shown in 
fig. 6 

y horizontal direction as shown in 
fig. 6 

z vertical direction 

? ratio of actual damping to critical 

damping 
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stiffening ring. (The spring rate of each spring is 300 Ib/ia 
(5. 254x10^ N/m). Each spring was stretched to 100 lb (4. 448x10^ N), 
i5 percent) 



-X 

Figure 4. - Stabilization system, a set of eight tangential springs at- 
tached to X“frame. (Springs 1, 2, 5, and 6 have a spring rate of 
200 Ib/in. (3.503x10^ N/m); springs 3, A, 7, and 8 have a spring rate 
of 100 Ib/in. (1. 751 x 104 N/m). Each spring was stretched to 100 lb 
(4w 448 x102 N), ±10 percent) 
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Figure s. - Typical mounting of piezoelectric-type accelerometers. 
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Figure 7. -Mounting of ioad ceii 
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Figure 9, - Electromagnetic shaker and X-frame. 
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(w) Frequency, 94.663 hertz. 
Figure 13. - Continued. 
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(b) Frequency, 47 hertz; damping ratio, 0.016. 



(c) Frequency, 72 hertz; damping ratio, 0.0145. 


Figure 15. - Typical transient data used for logarithmic decay determination of damping ratio. 
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